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ON THE USE OF THICK-AIRFOIL THEORY TO DESIGN AIRFOIL FAMILIES 
IN WHICH THICKNESS AND LIFT ARE VARIED INDEPENDENTLY 

By Raymond L. Barger 
Langley Research Center 

SUMMARY 

A method has been developed for designing families of airfoils in which the mem- 
bers of a family have the same basic type of pressure distribution but vary in thickness 
ratio or lift, or both. Thickness ratio and lift may be prescribed independently. The 
method which is based on the Theodorsen thick-airfoil theory permits moderate variations 
from the basic shape on which the family is based. 

INTRODUCTION 

In the design of airfoil sections, it is often desirable to vary the thickness or lift 
of a basic airfoil having favorable characteristics while maintaining as nearly as possi- 
ble the same type of pressure distribution. A rough approach (ref. 1, ch. 4) to this prob- 
lem is to design the mean line by thin airfoil theory and then superimpose on it the thick- 
ness distribution of a symmetric airfoil whose pressure distribution is known. 

Another approach is to start with an airfoil having desirable characteristics and 
scale the ordinates by a constant factor to obtain an airfoil with a different thickness. 
Airfoils designed in this manner have several shortcomings, as described in reference 2. 
Reference 2 then describes a different procedure which utilizes the thick-airfoil theory of 
Theodorsen. This procedure produces airfoils that differ from the basic shape both in 
thickness and in lift by essentially the same ratio. However, it is sometimes desirable 
to change the thickness of the airfoil without changing the lift, or to change both lift and 
thickness but in different proportions. The purpose of this paper is to describe a method 
for generating such variations within the context of the Theodorsen theory. 

SYMBOLS 


a arbitrary scale parameter 

C constant 
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c 


airfoil chord 


c jf airfoil lift coefficient 

Cp pressure coefficient 

R ,<p polar coordinates of exact- circle transformation of airfoil 

r ,6 polar coordinates of near-circle transformation of airfoil 

V undisturbed free -stream velocity 

v local velocity 

X = x - x Q 

x,y rectangular coordinates of airfoil in physical plane 

x 0 value of x at nose of airfoil 

a angle of attack 

e function relating angular coordinates of near-circle and exact-circle airfoil 

transformations 

value of e at airfoil nose 

e T value of e at airfoil tail 

6*,<p* dummy variables of integration 

t P 1 function relating radial coordinates of near-circle and exact-circle airfoil 

transformations 

arbitrary small constant 

average value of \f/ 

Prime superscript indicates derivative with respect to 8. 

2 

PRECEDING PAGE BLANK NOT FIUviED 



THEORETICAL CONSIDERATIONS 


The Theodorsen airfoil theory (ref. 3) involves a conformal transformation of the 
airfoil into a shape approximating a circle. (See fig. 1.) In the transformed plane this 
near circle is described by polar coordinates r ,9. The approximate circle is then 
transformed into an exact circle having coordinates R,(p. The function \p and the con- 
stant \p are defined by r = ae ^ and R = ae^o, respectively. Theodorsen shows that 
ip Q is the average value of 


and that the functions \p - ip Q and e = <p - 9 are related by the equations 


0) 




ip cot 


0 - 0 


2n Jq 


d <[>' 


0 - 0 ^* 

0 2 


e cot 


d0>’ 


( 2 ) 


( 3 ) 


The function 0/ is directly related to the airfoil coordinates in the physical plane by 


x = 2a cosh rp cos 9 


(4a) 


y = 2a sinh ip sin 9 


(4b) 


In practice i p is obtained as a function of 9 from equations (4), and then a func- 
tion e(6) is obtained from ip(Q) byreplacing <p* and <p in equation (2) with 6* and 
6, respectively. Henceforth in this paper e will refer to this function e(0) which is a 
close approximation to the exact function. The error involved in this approximation arises 
primarily as a tiny adjustment in e (but not in e’) in the velocity equation: 

v [sin(a + e + 8) + sin (a 1 + ej)](l + e’)e^° 

■jjf = . ■■ — - (5) 

y(sinh20/ + sin2(?)(l + ip's 


This discrepancy is normally very small compared with that arising from viscous effects, 
and so the approximation is considered to be adequate for the purposes of this analysis. 

Reference 3 also demonstrates that the angle of zero lift is -e T and the ideal angle 
of attack is -^j + £t)/ 2. These fundamental relations from Theodorsen’s analysis yield 
an interpretation of the airfoil family design procedure described in reference 2 which con- 
sists of simply multiplying the function by a constant C. According to equation (3), ip is^ 
thereby changed by a factor of C. Expanding equation (4b) in powers of ip yields 
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y = 2a sin 9 

Inasmuch as xp consists of \p Q plus small perturbations about i// Q , it is seen that the 
primary effect on the geometry of the airfoil is to change its thickness. 

It is seen that multiplying e by C automatically multiplies e T and e N by C. 
Consequently, both the angle of zero lift and the ideal angle of attack are multiplied by C 
and thus the design lift coefficient is changed by a factor of C. 

A second type of change is now considered: that of adding a constant ip a to xp. 

The constant \p Q is thereby changed by xp a but xp - ip Q is unaffected, as is its con- 
jugate function e. Therefore, the thickness of the airfoil will be changed but not its lift 
parameters. This possibility was mentioned in reference 4, 

Finally, consider the case where it is desirable to change both the thickness and the 
lift but in different proportions. This variation may be accomplished by combining the two 
procedures: first change the lift and thickness in the same ratio by multiplying i p and e 

by the same constant, and then add a constant to \p to adjust the thickness to the desired 
value. 



EXAMPLES 

In order to check the usefulness of the methods described in the previous section, 
several calculations were performed. The basic airfoil on which the family is based was 
very similar to that described in reference 5. The favorable characteristics of this air- 
foil are predicted by the method of reference 6 which includes a boundary -layer calcula- 
tion. The well-known effect of the boundary layer is to reduce the lift, as can be seen by 
the comparison in figure 2 of the viscous pressure distribution with the inviscid pressure 
distribution. In this case c £ was reduced from 0.455 to 0.381. 

It appears reasonable that although the predicted lift for the viscous case is signifi- 
cantly less than that for the inviscid case, a prescribed change in the inviscid lift will 
change the lift for the viscous case in approximately the same ratio. The assumption that 
is made here is that the boundary layer has approximately the same effect for all members 
of a family. Of couse, a slight error results from the assumption. The performance cal- 
culations in the following examples have been computed by the viscous theory in order to 
give a more realistic assessment of the usefulness of the design theory. 

Figure 3 shows the calculations for a family in which both the lift and the thickness 

were varied in the same ratio. The 17-percent-thick (Thickness = Q 1? \ airfoil . g th 
, . jit-,, \ Chord / 

basic model. For the 19-percent-thick airfoil the design calculation predicts a lift coef- 
ficient of 0.426 compared with the value of 0.429 obtained with the viscous calculation. 
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For the 15-percent-thick airfoil the design and viscous calculations yield 0.336 and 0.338, 
respectively. The difference amounts to less than 1 percent in both cases. The and 
e functions for this family are shown in figure 4. 

The airfoils shown in figure 5 vary in thickness from 15 percent to 19 percent but 
are all designed to have the same viscous lift coefficient 0.381. In this case the viscous 
calculation gives, for the 19-percent-thick airfoil, c^ = 0.388, and for the 15-percent- 
thick airfoil, Cf = 0.373, a discrepancy of about 2 percent in each case. These airfoils 
all have the same e function as the 17 -percent-thick airfoil. (See fig. 4(b).) The ^ 
functions for the 19- and 15-percent -thick airfoils are displaced by ±0.021, respectively, 
from that of the 17-percent-thick airfoil. 

Figure 6 shows another variation based on the same 17-percent-thick airfoil. In 
this case the thickness was increased by a factor of 1.059 to 18 percent, whereas the lift 
coefficient was designed to increase by 15 percent to 0.438. The lift coefficient actually 
obtained is 0.447, an increase of 17 percent. 

Naturally, it is not possible to design for a much larger increase in lift than thick- 
ness. Such an attempt results in an airfoil with zero thickness in the cusp region near 
the trailing edge. 


CONCLUDING REMARKS 

A method for designing families of airfoils has been described. These families are 
based on utilizing an initial airfoil with a desirable type of pressure distribution and then 
varying the thickness and/or lift while maintaining the same basic shape of the pressure 
distribution. The method is based on the Theodorsen transformation and so is not limited 
by the approximations of thin-airfoil theory. However, attempts to design large variations 
from the original airfoil tend to result in unsatisfactory shapes. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , February 14, 1974. 
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Figure 6.- Pressure distribution for an 18 -percent-thick airfoil design from the basic 17 -percent- thick airfoil 

with a design lift increase of 15 percent. 



